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We show that the Plebanski based approach to transformation optics overlooks some subtleties 
in the electrodynamics of moving dielectrics that restricts its applicability to a certain class of 
transformations. An alternative, completely covariant, approach is developed that is more generally 
applicable and provides a clearer picture of transformation optics. 
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The concept of transformation optics enables coordinate transformations on electromagnetic fields to be physically 
. realized as complex materials. This idea, applied to electromagnetic fields, originated with Pendry who pointed 
out the specific relationship between spatial transformations and material properties, and demonstrated how it could 
be used to create novel devices. Earlier, closely related work by Greenleaf et al. 12] developed a similar concept for 
Q '■ electric current flow and applied it to impedance tomography. The original descriptions of transformation optics were 
formulated for purely spatial transformations [3, 0] • Based on a theoretical foundation developed by Plebanski Q 
and De Felice Leonhardt and Philbin Q described transformation optics in terms of differential geometry and 
^ provided a recipe for the material constitutive relations for transformations that involve both space and time. 
^Lj, However, the Plebanski/De Felice based formulation is non-covariant, and part of the interpretation ascribed by 
' Leonhardt and Philbin overlooks some subtleties of electrodynamics in moving dielectrics, restricting its validity to a 
Cl^i certain class of transformations. These arise from a non-covariant matching of the constitutive relations in a curved, 
vacuum space-time with the constitutive relations for an isotropic moving dielectric. We show that such a matching 
\ cannot be made for a more general, anisotropic moving dielectric. It is possible to make such an identification with 

■ a stationary (i.e. at rest with respect to the frame in which the fields are defined) dielectric with non-vanishing 
OO magneto-electric coupling, but this can only be identified with a velocity if the material is isotropic. Moreover, the 
^-H , applicability of these non-covariant formulations to complex initial materials is limited. Although special cases have 

been addressed in the literature Q, a completely general formulation is lacking. 

Here we present a completely covariant approach that is valid for all transformations, materials, and material 
\^ ' motion. Based on a modern, completely covariant, description of classical electrodynamics, it allows for a more 
, rigorous and clear explanation of the concept behind, and interpretation of, transformation optics. This approach 

■ attempts to provide a systematic framework for the study of transformation optics and could be especially useful for 
applications involving moving objects, general relativistic corrections. 

The paper is organized as follows: In Sec. |lT] we review the completely covariant theory of electrodynamics in 
modern language. Section IIIII discusses electrodynamics in moving dielectrics and the differences between isotropic 
and anisotropic moving dielectrics; in particular how the non-covariant equations differ for isotropic and anisotropic 
materials. Section llVI outlines the Plebanski based approach, and explains the subtleties involved when taking into 
account the discussion of Sec. IIIII Section |V] describes the concept behind transformation optics and presents an 
interpretation consistent with the development of Sec. [Hi The main results for applications in transformation optics 
are Eqs. (|29p and (|30p. A variety of examples are given in Sec. I VII in order to both illustrate the details of the 
procedure, and to demonstrate agreement with known results obtained through other means. We conclude with Sec. 
IVIII and discuss possible extensions of this work to more general materials. 
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II. MODERN ELECTRODYNAMICS 



In modern parlance, classical electrodynamics is formulated in very geometric terms. In fact the homogeneous 
Maxwell equations are simply constraints on the fields imposed by the space-time geometry, and are a purely geometric 
condition. The inhomogeneous vacuum equations arise quite naturally from an action principle and a dual operation 
possessed by a manifold equipped with a metric. We follow the development and sign convention of Ref. js|, and make 
use of the Einstein summation convention, where repeated indices are summed over. The metric ga/s (or it's inverse 
g"^) lowers (raises) indices, and for the speed of light we use c = 1. For more details on the geometric foundations of 
electrodynamics the reader may consult a myriad of excellent sources such as p^| - [T2j . 

To make this more concrete, begin with the connection covector A — (the covariant version of the contravariant 
4- vector potential ^^), the field strength F = Fp^i, is an antisymmetric tensor defined as the exterior derivative (a 
generalization of the differential of a function) of A, 

F = dA, F^, = A,,^ - A^,y (1) 

where the comma indicates a derivative. The components of F can be represented as a matrix, that for Minkowski 
space-time in Cartesian coordinates have values 
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The field strength tensor F encodes some information about the fields, namely the electric field strength and the 
magnetic flux. On the other side of this coin, the magnetic field strength and the electric flux are encoded in another 
tensor G, called the excitation tensor. The components of G can also be represented as a matrix, that for Minkowski 
space-time in Cartesian coordinates have values 
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In much of the existing literature the matrix expression given for G is not Eq. ([3]) but rather has Ha and Da 
interchanged. However, this obscures the physical nature of the fields, where, for example, the electric flux must be 
the purely spatial part of G. The origin of the difference in definitions of G and further implications of this difference 
will be discussed elsewhere [l^. It is sufficient to know that, as the proper dual of F, G exists on an equal tensor 
footing and transforms as a tensor rather than as a "tensor density" , an unnecessary complication. 

The completely covariant form of Maxwell's equations are 

dF ^ 0. (4) 

and 

dG = J (5) 

where J = Ja/B-y is the covariant charge-current 3 index tensor. 

There remains one piece to the electrodynamics puzzle: an equation describing how F and G are related. The 
most natural way to define such a constitutive relation uses an operation called the Hodge dual, that on a four 
dimensional manifold with a metric g = takes a covariant tensor and returns another covariant tensor such that 

(*F)^. = IVWl^f.r^cpg^'^g^'F^S: (6) 

where tf^yap is the totally antisymmetric Levi-Civita symbol. This means that -k can be represented as a tensor with 
components 

V = ^v^^"^-p•9'^''5''^ (7) 
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and contains information about the space-time. The components of *F can also be represented as a matrix, that for 
Minkowski space-time in Cartesian coordinates have values 
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The constitutive equation is 



which in component form reads 
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G = x(*F), 



(8) 



(9) 



(10) 



The tensor x contains information on the dielectric material's properties including permittivity, permeability, and 
magneto-electric couplings, and can be thought of as representing an averaging over all the material contributions to 
an action that describes a more fundamental quantum field theory. To retain the desired properties and usual notions 
of G and F, x must be independently antisymmetric on its first two and last two indices, and in vacuum x(*F) = ★F. 
Additional symmetry conditions are frequently impose d IT3. fisl for various reasons not considered here because we 
would like to be able to describe very general materials |l6l.ll7l|. but the last condition is sufhcient to uniquely specify 
all components of % for the vacuum, they are 



iXva 




(11) 



Equation expresses % as a matrix of matrices, the first two indices of Xap^'^ §i^^ component of the large 

matrix, which is itself a matrix described by the second set of indices. The component values of x for the vacuum 
are unique and independent of coordinate system. For a more general material, the component values can easily 
be determined by simply matching the results of the constitutive equation G — x(*F) with the usual flat-space 
constitutive relations in a particular coordinate system. The components of the constitutive equation provide a set 
of six independent equations that can locally be collected in the form 



H = fr^B + ii*E, D = e*E + i2*B. 
where we use the notation a to denote a 3 x 3 matrix. Rearranging these to 

B = fiH + iiE, D = eE + 73/? 



(12) 



(13) 



gives the more familiar representation for the constitutive relations. These three-dimensional representations of the 
completely covariant Eq. Q are essentially equivalent, and it is a simple matter to switch between them using the 
relations 



72 M7i 



71 



-A*7i 



72 = 72 M- 



(14) 



However, one should be aware that these 3x3 matrices are not tensors but simply components of x that have been 
collected into matrices. To recover the usual component relations of Eq. (IT^ the constitutive equation G = x(*F) 
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allows us to make the identification 
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where the * indicates entries that are antisymmetric on either the first or second set of indices on x. 



(15) 



7^ 



III. ELECTRODYNAMICS OF MOVING DIELECTRICS 

Consider a dielectric material moving with uniform velocity with respect to the lab frame. The material properties 
of the dielectric are known in the material frame ~ the frame at rest with the material - and are described by 
permittivity e, permeability jl and magneto-electric couplings 71 and 72. For a field tensor F measured in the lab 
frame, the components of F relative to the material frame are found by changing to the material frame with a Lorcntz 
transformation, L ,9.. ,18.] 



F' = L(F) ^ F^.p. 

Inside the material, G' — L(G) is now related to F' by the usual constitutive equation G' = x(*F'), or 

L(G) = x(L(*F)), 



(16) 



(17) 



resulting in expressions involving combinations of D and H in terms of combinations of E and B. On the other hand, 
the inverse transformation of G', 



G^L-Mx(L(*F))] 



(18) 



gives direct expressions for the components forming D and H . Furthermore, since the constitutive map G = x(*F) 
is totally covariant, a third representation can be obtained by transforming \ from the material frame back to the 
lab frame 



G=(L-i(x))(*F). 



(19) 



Equations p7p - (|19p provide equivalent representations for the constitutive relations in a moving dielectric, are 
the most convenient notation for dealing with electrodynamics in a 4-dimensional, covariant way, and are valid for 
arbitrary space-times and transformations. If one insists on expressing these four dimensional relations in terms 
of 3-vectors, then for the special case of a dielectric material moving through the Minkowski vacuum with uniform 
velocity. 



H+ (D X f3 



B+ (E X 13 



D 



H X (3 



B X 13 



71 



■ 72 



B X f3 



B+ (E X /3 



(20a) 
(20b) 



These equations are typically shown in the literature only for the restricted case where the material is isotropic, 
e = e, fi = ^, and has no magneto-electric coupling, 71 = 72 = (see, e.g. Ref. [l3l)- For this restricted case, the low 
velocity limit is 



D = sE+{efi-l){l3 X H), 



(21a) 
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B = fiH - {e^i-l)il3 X E). (21b) 

While these isotropic, low-velocity results may be the most familiar, they are not correct when e and are matrix 
valued. Indeed, following the same procedure leading from Eqs. (|20|) to Eq. pT|) with matrix- valued e and and 
with — j2 — 0, leads instead to 



D = eE + e{(3 X jlH] - [13 X H] , (22a) 



B = fiH + fi (^eE X 13^ - (^E X . (22b) 

Matrix multiplication is not commutative with the cross product. The difference in going from an isotropic material 
to a anisotropic material may seem trivial at first - and this is true when one maintains covariance, as in Eq. (|20|) - 
but by giving up covariance in passing to the low-velocity limit the results change significantly. Even diagonality of 
e and fi~^ is not sufhcient, one may readily show that recovering Eqs. (|2ip requires pure isotropy. 



IV. PLEBANSKI BASED TRANSFORMATION OPTICS 

Plebanski Q previously studied the propagation of electromagnetic waves in gravitational fields, and much of the 
current analysis in transformation optics is based on his results . Consider the propagation of electromagnetic waves 
in a vacuum space-time described by an arbitrary metric (/a/j. Using g^iB and % = Xvac in Eq- (IH), the constitutive 
relations for this vacuum space-time can be written 

Da = -^g'^'E, -f eabc—H,, (23a) 
9oo 9oo 



Ba = -^g''''Hb-eabc — Ec. (23b) 
500 goo 

Written this way these equations are not, as cautioned by Plebanski, covariant. Notice that while they still sum 
over repeated indices they do not conserve index type. Comparing with the results for a moving dielectric in the 
low velocity limit, Eqs. (|21|) . it would appear that a dielectric moving with low velocity through vacuum Minkowski 
space-time with permeability and permittivity 



and (scaled) velocity 



(24) 

goo 



gob fryr\ 

Vb = — , (25) 
500 



will reproduce the same field relations as the arbitrary vacuum space-time. 

However, as we have seen in Sec. Illli Eqs. (PT|) are only valid when the material is isotropic, therefore this approach 
is not valid in general. The problem stems from attempting to match the covariant Eqs. (|23l) with the non-covariant 
Eqs. (|2ip . It is clear that such a matching cannot be made with the more general Eqs. (|22|) . This method, therefore, 
only works when the material is isotropic (in which case Eqs. (|22l) reduces to Eq. (|2ip ). or when gob = for all b (in 
which case the velocity must be assumed zero). 

Despite these restrictions, this approach frequently gives, or appears to give, correct results. This is because there 
is a case where the covariant Eqs. can be matched with a covariant expression in dielectrics. Remember that 
the general constitutive relations for a stationary (at rest with respect to the frame in which the fields are defined) 
dielectric can be written in the form of either Eqs. (IT2]) or Eqs. (IT3)) . which are related by Eqs. ([T4|) . Matching Eqs. 
([T3| with Eqs. ([23| . one could identify e and jl as in Eq. ([24l) but with magneto-electric coupling 



72 =7i = -eabc — . (26) 
500 

It is important to realize that when trying to match Eqs. (|23l) with Eq. (|13p . that the E, B, etc. found in Eq. 
are not actually vectors but components of either 1-forms or 2-forms, and that interchanging a 1-form with a vector is 
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only trivial in Minkowski space-time. It is also important to realize that this magneto-electric coupling cannot always 
be simply interpreted as a velocity. Indeed, returning to the more general Eqs. (|20p one can find the low velocity 
limits (again in Minkowski space-time) 

D^eE + i2H - e (fiH x /3) + (i? x /?) - e (71^ x /?) -|- 72 [sE x /s) + 7-3 (fai? x , (27a) 

B = ilH + iiE + jj [eE xp^-[Exp^+jl {-j^H x /s) - 7'i x /s) - 71 (71^ x , (27b) 

from which it should be clear that the only way to identify a magneto-electric coupling with a velocity is to replace a 
stationary magneto-electric material with an isotropic moving dielectric with no magneto-electric coupling. Instead, 
it is better to employ a fully covariant procedure based on the theory as developed thus far, which is not only valid 
for a very general class of materials and arbitrary dielectric motions [s!] , but more clearly distinguishes between the 
material and space-time contributions, and is also valid for arbitrary background space-times, such as the weakly 
curved space-time around Earth. 



V. TRANSFORMATION OPTICS 



Due to the apparent limitations and non-covariance of the Plebanski/De Felice based approach to transformation 
optics, we ask whether there is a more general approach. By following the covariant development of electrodynamics 
in the previous sections we find that there is. In an attempt to more clearly tie the concept to the method, consider 
the following interpretation of transformation optics. 

Start with an initial space-time manifold (M, g,*), field configuration (F, G, J), and material distribution x, where 
dF = 0, dG = J, and G = x(*F). Imagine now a map T : M ^ M C M that maps M to some image M and 
transforms the electromagnetic fields in some smooth way to a new configuration (F, G, J). Because the underlying 
space-time is physically unaltered the manifold is still described by (M, g,*). But for the new field configuration to 
be physically supported there must exist a new material distribution %. Thus on the image M we must have dF — 0, 
dG = J, and G = x(*F). Such a transformation could, for example, map M to an image M that contains a hole, 
i.e. a region from which the fields will be excluded, as in the case of an electromagnetic cloak [H, [l^. 

Using the inverse T of the map T (but see the remarks below), we can relate the initial F and G to the final F 
and G by an operation called the pullback of T, which we denote as T* ■ Schematically, this means that we have the 
relation 



G = r*(G) = r*(x(*F)) = x{*T*{F)). 
This can be solved for x as a function of x e M to find [l^ 



(28) 



T{x) 



r{x) 



T7] 



(29) 



In Eq. A is the Jacobian matrix of T, A^^ is the matrix inverse of A, both A and A^-'^ are evaluated at x, and in 
solving for x wc have made use of the fact that on a 4-dimcnsional Lorcntzian manifold, acting twice with -k returns 
the negative, i^-kF ~ — F. 

Note that the initial material parameters must be evaluated at T{x), but if the initial space-time is vacuum, then 
since Xvac* = *, 



T{x) 



(A-i)^(A-i)^ 



(30) 



where one * is evaluated at T{x), and everything else is evaluated at x. The prescriptions of Eqs. ([^5]) or ([50)) are 
meaningful only for points x M . So for transformations such as the electromagnetic cloak, where there is a hole 
in M, the material parameters inside the hole are unspecified and completely arbitrary. In this way, any uncharged 
material may be hidden inside the cloak without affecting the behavior of the fields outside. 

Equation ([^5]) is the core of transformation optics. Start with a given space-time with metric g and associated dual 
and with given dielectric material properties described by the tensor x- The initial space-time may be Minkowski 
and the initial dielectric may be the vacuum, but this is not necessary. We imagine a transformation that changes 
the fields in some way. We ask what x is required to physically achieve such a transformation. The answer is given 
by Eq. (EH). 
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Two final remarks are in order. One is to note that the crucial map is not T, but rather T- So one could just as 
well start with some map that here we label T, and completely ignore T . This is important because a given T may 
not have an inverse, so we must assume that the given T is well defined on its domain, M. But for practical purposes 
T can be thought of as the inverse of T. Secondly, while Eqs. (23) and (15(11) would be quite difficult to evaluate by 
hand, a modern computer algebra package makes the calculations almost trivial. 

VI. EXAMPLES 

The completely covariant approach to transformation optics developed above provides a concrete and powerful 
framework for analysing any desired configuration of fields and linear dielectric materials in any space-time and with 
any relative velocities. This section presents some examples to illustrate the usefulness of the covariant approach 
described above and show that it recovers previous results obtained through other means. In these examples we use 
a notation whereby the point (t' ,x' ,y' , z') is a point in the original, vacuum, space-time, while a point {t,x,y,z) is 
a point in the space-time plus material. We seek to describe the material properties as a function of the unprimed 
coordinates, x, y, z). 

A. Spatial Inversion 

For a simple inversion of one spatial coordinate, the transformation 

Tit', x', y', z') = it, X, y, z) - {t' , -ax', y' , z') (31) 
forms the basis for a "superlens" [13, [HI and has the inverse map 

Tit, X, y, z) = it', x', y', z') ^(^t,-^,y,zy (32) 

The Jacobian matrix of T is A"^ = diagil, — 1/a, 1, 1). Using this in Eq. (jSOj results in a x, that, upon comparison 
with Eq. ([T5|). leads to 

/-a \ 

(33) 




B. Temporal Inversion 

On the other hand, inverting the time coordinate with the transformation 

Tit', x', y', z') = it, X, y, z) = (-at', x' , y' , z') (34) 

leads to 

/-a \ 

e = A= -a . (35) 
\ -a) 

Of course, replacing a — ^ —a in the previous example corresponds to a simple stretching of the time coordinate, and 
the result is the same as Eq. (1351) with a — J> —a. Both this example and the previous highlight the ease with which 
the completely covariant method handles coordinate inversions - a concern raised in Ref. |21|] . 

C. Spatial Dependent Time Transformation 

Next consider a transformation that has also been analysed in more detail from the perspective of creating frequency- 
altering linear materials, using other methods [2^ . namely 

Tit', x', y', z') = it, X, y, z) = ((ax' -t- h)t' , x' , y' , z') (36) 
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In this case the inverse map is 



T{t,x,y,z) = (t ,x ,y ,z ) 



(37) 



and the Jacobian matrix of T is 



/ ax+b (ax+by^ ^ ^\ 

10 

10 

V ly 



(38) 



The next step is to calculate x using Eqs. ([3D]). The results can be extracted from identifications with Eq. ([TS]), and 
written in the representation of Eq. as 



fax + b 

(^(ax+b}^+at)(^{ax+b)^-at) 







V 















(^(ax+b}^+at)(^{ax+b)^-at) 
(ax+bj^ 



1 



M = 



ax + b 



'1 0> 
1 
,0 1, 



(39) 



71 = -72 




{ax + 6)2 



When rewritten in the more common representation of Eq. (|13p they become 



'1 ON 

e = li= (ax + b) \0 1 
.0 1, 



at 



7i = 72 



ax ■ 



^0 
0-1 
,0 1 



(40) 



This example highlights the slight difference in the meaning of what we are calling e*, p,~^, 71*, and 72* and the 
interpretation of these quantities based on the representation of Eq. (|13p . This difference arises only when non-zero 
magneto-electric coupling terms are present, and it is easy to switch between the two representations using Eq. (jl4p . 
While the computation and manipulation of large matrices like x rnay see rather daunting, the use of a modern 
computer algebra package makes the calculations almost trivial. 



D. Time Dependent Spatial Transformation 



Having looked at the behaviour exhibited by a spatial dependent time transformation, it is natural to next enquire 
about a time dependent spatial transformation. Let the transformation be 



The associated map T is 



and the Jacobian matrix of T is 



T(t', x', y\ z') = (t, X, y, z) = (t', {at' + b)x' , y' , z') 



at + b'^''}' 



T{t,x,y,z) = {t\x',y',z') ^ t 



(I 0\ 

ax 1 n (-) 

A/i — {at+b)'^ at+b " " 

10' 

\ 1/ 

Proceeding as before and expressing the results in the representation of Eq. ([T3| , 

0' 



e = fi 




{at + 6)4 - a'^x 





{at + b)^-a^x^ [I 5 ~o 



ax{at + b) 



(41) 



(42) 



(43) 



(44) 
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It is interesting to compare the last few examples. In the cases of a simple scaling of t versus a spatially dependent 
transformation of t, the resultant materias have similar characteristics in that e is isotropic and scales like t, while 
the mixing of space and time introduces a magneto-electric coupling. According to the discussion in Sec. IIV| because 
the resulting material is isotropic, this magneto-electric coupling, if small enough, could arise from a small material 
velocity. So the resulting material could be considered as being cither at rest and possessing a magneto-electric 
coupling, in motion with no magneto-electric coupling, or some combination of the two. 

Turning to the cases of a simple scaling of x versus a time dependent transformation of x, Sxx still scales like x, but 
Eyy = Ezz have totally different behaviours in each case. Again the mixing of space and time introduces a magneto- 
electric coupling, but this time, because the resulting material is not isotropic, these magneto-electric couplings cannot 
be interpreted as a material velocity. 

To verify the material parameters obtained here, a validation calculation must be performed. Such a calculation 
should be based on a combined space-time transformation in a complex initial material to demonstrate the full 
generality of our approach. Such a calculation, however, requires an in-depth analysis for both the initial medium 
and the transformed medium that is beyond the scope of this paper, which is intentionally a focused description and 
formulation of the new idea. However, such work is in progress. 



E. Square Cloak 



To demonstrate agreement with previous results, consider now the well-known, but non-trivial, example of a square 
cloak, studied by Rahm et. al. fl^. In this example the transformation is a map T piecewise defined by 



T{t\ x\ y', z') = {t, X, y,z)= [ t, x' 



S2 



S2 



si,y 



S2 



S2 



fl 

x' 



(45) 



for (0 < a;' < S2), (— S2 < y' < S2), \y'\ < \x'\, and \z'\ < 00, and T{t' ,x' ,y' , z') = {t,x,y,z) elsewhere. Again, we are 
really only interested in the inverse of this map, which we take to be piecewise defined by 



T{t, X, y, z) = {t', x', y', z')= it 



S2 , s S2{x-sx)y 

S2-S1 {S2-Si)X 



(46) 



for (si <x< S2), {—S2 <y < S2) with \y\ < \x\ and \z\ < 00, and T{t, x,y, z) — [t' ,x' ,y' , z') for x > 82- This map is 
not defined for a; < si, so x is undetermined and therefore arbitrary in this region. The Jacobian matrix of 7" is 



A" — 



'I 



Vo 





S2 



ySlS2 






S2(X — si) 



' 





(47) 



^(S2-Sl) X(s2-Si) 

01/ 



Turning the crank on Eq. (|30p . and comparing with Eq. (jlSp the components of e and p. can be identified and gathered 
in matrix form as 



M I x^{x — s-i) 








jlix-si) 

X{si-S2y-' / 



(48) 



and 71 = 72 = 0, which are exactly the results obtained in Ref. [l^. Note that the domain of 7", {si < x < S2), 
(—§2 < y < 82) with \y\ < \x\ and \z\ < 00, dictates also the domain of e = /i, so the material distribution is defined 
piecewise, with this non-vacuum piece being located exactly where it is desired. 



VII. CONCLUSIONS 



Having shown that there may be important limitations to the approach to transformation optics based on the 
constitutive equations originally presented by Plebanski [H, we ask whether it is possible to construct a completely 
covariant and more general approach to transformation optics in linear materials. The main result of this construction, 
Eqs. and (|5D|) . represent the core of transformation optics. 

While the Plebanski based approach remains useful and valid for a wide class of transformations, the benefit of this 
approach is that it is valid for arbitrary background space-times, arbitrary initial material, arbitrary material motion. 
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and for arbitrary transformations - including those that mix space and time as. At first sight the totally covariant 
approach seems to have the drawback of being more difficult to work with since it involves transforming a matrix 
with 256 components. But with the aid of a computer algebra package these kinds of computations become trivial, 
and the recipe given by Eqs. (23) or ([50]) is straightforward and applicable in all situations. 

Several examples have been discussed that illustrate the material parameter features that might be expected from 
various types of transformations. Particularly interesting are the examples of transformations that mix space and 
time. It is found that a generic feature of these types of transformations is the appearance of magneto-electric 
coupling terms. For transformations that have only a spatially dependent time transformation, the resulting material 
is isotropic, and the magneto-electric coupling could be reinterpreted as a velocity, in accordance with the Plebanski 
based method. However, if the transformation mixes a spatially dependent time function with an additional spatial 
function of any type, or if the transformation is of a time dependent spatial function type, then the resulting material 
is not isotropic, and the magneto-electric couplings cannot be reinterpreted as material motion. 
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